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Diffusion and reaction of initially separated ions and in the presence of counter ions A'^ 
and B~ is studied. The dynamics is described in terms of reaction-diffusion equations obeying local 
electroneutrality, and the time-evolution of ion-concentrations is determined. We find that, in the 
absence of reactions, unequal mobility of ions generate nontrivial features in the macroscopically 
observable concentration profiles. Switching on the reaction A~ + B'^ — > C leads to the formation 
of a localized, diffusive reaction front and one finds that the properties of the front (e.g. effective 
diffusion constant) are affected by the background ions. The consequences of this effect on the 
formation of Liesegang patterns is discussed. 

PACS numbers; 0.5.70.Ln, 45.70.Qj, 66.10.-x, 82.45.+Z 



I. INTRODUCTION 

The reaction-diffusion process A+B C has been dis- 
cussed for a long time. This conceptually simple process 
displays a rich variety of phenomena (nonclassical reac- 
tion kinetics Jp|,H , clustering and segregation |^,^, front 
formation [plra]j and, depending on the interpretation of 
A and B (particles, quasi-particles, topological defects, 
chemical reagents, etc.), it provides a model for a number 
of phenomena in physics, chemistry, and biology. 

In many cases of interest, A and B are ions {A~ and 
B^) and these ions are initially separated from each 
other. An example we shall discuss below is the forma- 
tion of Liesegang bands where an electrolyte A'A'^ 
diffuses into a gel column containing another electrolyte 
B^B^. The concentration of As is taken to be much 
larger than those of the B-s, thus the reaction front 
A^ + 5+ C moves along the column. An appropriate 
choice of reagents then leads to quasiperiodic precipita- 
tion (C —> D) in the wake of the front (Fig.l). 
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FIG. 1. Schematic picture of Liesegang phenomena. The 
correspondence with the notation in the text is given by 
A~ = Cr, A+ = H+ (outer electrolyte); B+ = Ag+ , 
B^ — NO^ (inner electrolyte); and D = AgCl (precipitate). 
The initial interface between electrolytes is at a:: = 0. The 
precipitation bands (shaded regions) emerge in the wake of 
the moving reaction-diffusion front (dashed line at Xf). 



In general, the background ions (A+ and B~) are ex- 
pected to play a role in a process described above. Nev- 
ertheless, the usual approach is to neglect them and con- 



sider only a contact interaction between neutral reagents 
A and B. This approximation is based on the argument 
that the background ions provide only screening and, fur- 
thermore, the screening length is much smaller than the 
scale of concentration variations relevant in the forma- 
tion of a macroscopic pattern. Although the argument 
sounds compelling, one should note that the background 
ions may generate macroscopic effects even if the screen- 
ing length is negligible. Indeed, if the mobility of one of 
the background ions (A+ in the Liesegang case) is much 
smaller than the other mobilities then the motion and 
the properties of the reaction front are altered. Since the 
properties of the reaction front are crucial in determin- 
ing the pattern ||9|-[ll| one expects that the presence of 
background ions gives rise to macroscopic changes in the 
observed patterns. 

Our aim with this work is to verify the above expec- 
tation and to investigate how the diffusion and front for- 
mation are affected by unequal mobilities of background 
ions. More precisely, we shall study the time evolution 
of ion-concentrations in the process 
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(1) 



where the reaction product C = A B is assumed 
to vanish from the system. The process starts at 
t = form an initial condition where the electrolytes 
A^A'^ and B'^I3~ are separated and their concentra- 
tions (a+, a^, 6^, 6+) are constant in the left {x < 0) and 
right [x > 0) half-spaces, respectively 



a~ (x, t ^ 0) = a+{x,t = 0) = aQe{-x) 
h+{x,t = 0) = b-{x,t = 0) = baOix) 



(2) 



where 9{x) is the step function. Such an initial state 
with ao ^ 6o is actually used in Liesegang experiments, 
and this choice is also motivated by the fact that in- 
vestigations of front formation from such initial state 
have proved to be instrumental in understanding the 
A + B ^ C process [||. 
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The study of motion of ions is not an easy task and we 
must simphfy the problem to make it tractable. We be- 
lieve, however, that our approximations listed below are 
appropriate at least for the description of the Liesegang 
experiments. 

1. It is assumed that the phenomena can be described 
by reaction-diffusion equations. This appears to be 
a correct assumption for reactions taking place in 
a gel where convection is absent. 

2. The screening length is assumed to be negligi- 
ble and screening is taken into account by enforc- 
ing local electroneutrality. At characteristic ion- 
concentrations — IM) present in Liesegang 
experiments, the screening length is indeed small 
(~ 10~^m) compared both to the characteristic dif- 
fusion length 10 
reaction zone (~ 10 
be found in Sec.|l[ 



to) and to the width of the 
^m). Further discussion can 



3. The concentration profiles are assumed to depend 
only on one spatial coordinate (x in Fig.l). Al- 
though a one dimensional geometry can be set up in 
experiments on Liesegang phenomena (the length 
of the gel column can be made much larger than 
its width), one should note that the finite extent of 
the sample in the transverse direction poses non- 
trivial problems with edge effects. It appears, how- 
ever, that these effects can be neglected since the 
final pattern is usually one dimensional to a good 
accuracy. 

4. The mobility of the reagents and of the background 
ions are, in general, different. For simplicity, we 
shall consider the case with one of the background 
ions having a significantly distinct diffusion coeffi- 
cient 



Dr, - Dh 



D^Da = D. 



(3) 



This is just a technical assumption to keep the num- 
ber of parameters small and, this also appears to be 
the most interesting case for Liesegang phenomena 
where bo- 



properties of the reaction front are calculated. Finally, 
the implications for understanding the Liesegang phe- 
nomena is discussed in Sec.^ 



II. EQUATIONS IN THE 
ELECTRONEUTRALITY APPROXIMATION 

In a medium such as a gel, the ions move by diffusion 
and, in the presence of an electric field E — — V(/3, the 
flux of ions ji is given by the Nernst-Planck relation |12] 



= •?i,diff + ■?i,drift = y^"^^ + —'^'^^ ) ■ (4) 

Here is the concentration of i-th ions of integer charge 
Zi, Di is their diffusion coefficient, and ifQ = RT/F is 
a constant combined of the temperature T, the gas con- 
stant i?, and the Faraday number, F . The potential Lp is 
determined from the Poisson equation 



A(/5 = > Zim 



(5) 



where eo is the permittivity of free space while £r is the 
dielectric constant of the system. 

An important quantity in ionic diffusion is the Debye 
length ty) that gives the characteristic length-scale asso- 
ciated with charge inbalances 



F^m 



(6) 



where uq is the characteristic scale of ionic concentra- 
tions. In a Liesegang experiment, one usually has no ~ 
10^'^M — IM and the process takes place in an aqueous 



medium [sr ~ 80). Thus ro 



10 



-10^ 



10 m and 



one can see that rjj is much smaller than the scale of 
the macroscopic pattern (e.g. the width of the bands, 
« 10~^m — 10~^to). As a consequence, one can use the 
electroneutrality approximation that consists of replacing 
eq.(||) by the constraint 



^2 ^i'^i = ■ 



(7) 



Once the above approximations are made one arrives at 
a problem that can be studied numerically and, in some 
limits, analytically. The process is now simple enough so 
that the numerical analysis is not hindered by computer 
time and memory problems, or by difhcultics arising from 
discretization. 

In order to arrive at the results, we shall proceed as fol- 
lows. First we discuss how to take into account the elec- 
troneutrality constraint in the reaction-diffusion equa- 
tions (Sec.|^. Then the case without reaction is stud- 
ied and we show that interesting concentration profiles 
emerge even in the pure diffusion process (Sec. III). The 
effects of reactions are considered in (Sec.tV) where the 



Denoting now the rate of reaction of the i-ih. ion with 
the others by Ri{{n\) and assuming that the reaction 
does not violate electroneutrahty i^iZiRi = 0), eq.(Q) 
together with the constraint (0) yields the following equa- 
tions for the time-evolution of the concentration fields 



dfjii = D. 



W-{n,£)] ~R^{{n}) 



(8) 



where the appropriately scaled electric field that arises 
from the electroneutrality constraint is given by 



J2i ZjDiVni 



(9) 
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It should be noted that there is an extra term in £" if a 
steady global current flows through the system. Such a 
current is not present in the Liesegang problem and, try- 
ing to keep the discussion as simple as possible, we shall 
assume that the global current is zero. 

For the process of actual interest (|l|), reaction takes 
place only between the ions A" and and their rate of 
reaction is given by ka~b'^ where k is the rate constant. 
Thus the above equations in a one dimensional geometry 
take the form 



with 



dta- = D 


dla 




- ka-b+ 


(10) 


dtb+ = D 


'dlb+ 


~d^{b+£) 


- ka-b+ 


(11) 


dta+ = D 


'dla+ 


-d^{a+£) 




(12) 


dtb- = D 


'dfb- 


+ d^{b-£) 




(13) 


^ _ DdA- 


-a" - 


^b+-b-)^ 


- Dd^a+ 


(14) 


D{a- + 


b+ + b-) + Da+ 



Equations (|l0|-p^ together with the initial conditions (||) 
provide the mathematical formulation of our problem. 

Before turning to the solution of the above equations 
let us mention that the diffusion-reaction problem of ions 
in one-dimensional geometry can be tackled numerically 
without assuming the electroneutrality condition. The 
only difficulty is that the discretization of space must 
be on a finer scale than the Debye length and so, in 
the range of physical parameters where is exceedingly 
small, the calculation becomes impractical. One expects 
(and we verified it for some cases) that the solution of the 
full problem approaches the solution of the corresponding 
"electroneutral" problem as is decreased. 



CONCENTRATION PROFILES WITH NO 
REACTIONS 



III. 



Let us begin the analysis of eqs.(|lC|-|l4|) by considering 
the case of no reactions {k — 0) and let us further re- 
strict our study to the case of oq ^ 6o corresponding to 
the Liesegang initial conditions. The limit of bo = is es- 
pecially simple and treated in textbooks. In this case, the 
two ions A~ and A'^ must move together thus an electric 
field is generated that slows down the more mobile ions 
and accelerates the slower ions. The result is an effective 
diffusion with a diffusion coefficient DcS = 2DD/{D + D) 
@- 

The presence of a small amount of B-s (oq 3> bo) does 
not significantly change the motion of As. The ions A~ 
and A+ can separate now but only by an small amount 
that is compensated by the motion of J5-s. On Figj^, 
we can see the results for the case of bo/ao = 0.01 and 
D/D = 0.1 (slow background ions A+). 
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FIG. 2. Concentration profiles of the A and B ions [(2a) 
and (2c), respectively) and the electric field (2b) generated 
by them. The concentration is measured in O.IM while the 
unit of the electric field is V/m. The results are for the 
case of bo = O.Oloo with the diffusion coefficients given by 
D — lQ~®m^/s and D = Q.ID. Inset shows that region where 
relative separation of the A ions is significant. 



The electric field is mainly generated by the motion of 
the majority A ions and, in turn, this field is the deter- 
mining factor in the motion of the minority B ions. Since 
this field, shown in Fig.||b, moves the 4+ (A^) ions in 
the +x {—x) direction, similar effect is felt by the B ions. 
Indeed, as one can see in Figj^c, the B^-s are repelled 
from the region the A ions moved into, while the ions 
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I3~ are pulled through this region. As a result a region 
emerges where the ionic and diffusive drift of the B-s are 
in opposite directions. 

It should be noted that the profiles shown in Fig. 2 keep 
their shape in time. The pictures at t' are obtained from 
those at t = 1ft, by rescaling the x-axis by a factor \fV~fi- 
This numerical observation is the consequence of the fact 
that the Debye length is zero and the initial conditions 
(|) do not contain any length-scales. As a consequence, 
all the length-scales are diffusive lengths proportional to 
\ft. The above argument can be seen to work explictely 
in the limit Z) = where an analytic calculation gives 
the concentration profiles that can be expressed through 
error functions of argument xj \Jt. 

The complexity of concentration profiles shown in 
Fig. 2a and 2c suggest that if reactions are switched on 
between the ions A~ and i?+ then the emerging reaction 
front may be rather different from the case of neutral 
reagents. This is what we shall study in Sec.|^. 



above picture is modified as a result of the ionic charac- 
ter of the reagents. 

Eqs.( |lO[jr^ ) with initial condition (||) can be studied by 
straightforward numerical methods and one finds that 
the localized-diffusive-front picture does hold and, fur- 
thermore, the scaling properties (p^p^ also remain valid 
when the ionic interactions are switched on. The actual 
values of the parameters D f and cq, however, are affected 
by the presence of background ions. 

In order to understand how these results arise, let 
us begin with the numerical observation that the re- 
action front remains narrow even if the ionic interac- 
tions are switched on. Indeed, for characteristic val- 
ues of ao ~ IOO60 ~ IM, Da ^ Db ^ 10-i°m2/s and 
k w 10^° (Ms) ~i we find that the width is in the 
mesoscopic range 10~^m) at all times available in a 
Liesegang experiment. Thus, on diffusive length-scales, 
the reaction zone can be treated as a point (as in the neu- 
tral case) and one arrives at equations with no reaction 
terms 



IV. REACTION FRONT 



(18) 



The full reaction-diffusion process is described by 
eqs.(pii|-p^ and the solution of these equations with ini- 
tial condition (H) provides the description of the reaction 
front. Indeed, once the concentration profiles are known, 
the location and the time evolution of the production of 
A~ + _B+ — > C particles is given by 



R{x,t) = ka {x , t)b'^ {x , t) . 



(15) 



The properties of R{x, t) are well known for the case of 
neutral reagents (5 = 0) In that case, the reaction 
takes place in a narrow, moving region whose width is 
much smaller than the diffusive scales. The motion of 
the reaction zone is 'diffusive' characterized by a diffu- 
sion constant Df 



Xf = y/2Drt. 



(16) 



Another important feature of the front is that it leaves 
behind a density of C-s pOl] 



Co 



R{x,t)dt, 



(17) 



that is independent of x. 

The parameters D j and cq can easily be determined for 
the neutral case by exploiting the smallness of the width 
of the reaction zone. The reaction zone is replaced by a 
point where the diffusion equations are supplemented by 
boundary conditions and as a result the parameters Dj 
and Co are given as functions of Oq, bo, Da and Df, []5|,p^. 

The presence of a localized, diffusive front is an essen- 
tial ingredient in the theories of Liesegang phenomena 
[p|-|rH, and the parameters of the front (especially Df 
and Co) are known to influence the properties of the pat- 
terns. Thus the next step is now to find out how the 



The reactions are taken into account by the following 
boundary conditions at the front 



a (xf) = b+{xf) = 0, 

lia-(2:/)| = \jb+{xf)\ 



(19) 



The meaning of the above conditions is that the concen- 
trations of the reagents are zero at the front and the flux 
of ions and to the reaction zone are equal. 

Let us now suppose that Xf(t) and ni{x, t) are the so- 
lutions of the above equations (|l8-19) with the initial 
condition (||). Then one can easily verify that the front- 
position Xxf{t/\^) and the concentrations ni{x/X,t/\^) 
also solve the same problem for an arbitrary A > (note 
that the initial conditions do not contain any length- 
scale). Thus the functions ni(x,t) and Xf{t) must sat- 
isfy the conditions ni{x,t) — ni[x/\,t/\^) and Xf{t) = 
Xxf{t/\^). As a consequence, we flnd that the concetra- 
tion proflles obey the following scaling form 



ni{x, t) = <f>i 



(20) 



and the front moves diffusively even if the ionic interac- 
tions are taken into account 



Vt. 



(21) 



The above relations hip ([2l| ) defines the diffusion constant 
Df through Xf — y/2Dft. 

The scaling of the concentrations (^) together with 
equation (ph imply scaling for the electric field 



£{x,t) 



Vi \VtJ 



(22) 
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These scaling results allow us to investigate the pro- 
duction of the C particles. The number of the C-s arising 
in the reaction zone in a unit time is given by the flux 
of one of the reagents (e.g. j^-) entering the front. Ac- 
cording to eqs.(|2C-22) j^- at Xf is proportional to \/\/t 
and the velocity of the front decays in time in the same 
way, Xf ^ It follows then that the density of the 

C-s emerging in the wake of the front is a constant 



J a 
Xf 



const. 



Co 



(23) 



The results (|20(-^3|) given by the above analytical ar- 
gument have been confirmed by computer simulations. 
An example of such numerical calculation can be seen in 
Fig.|. 
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FIG. 3. Concentration profiles of tlie A and B ions wlien 
the reaction is switched on. The initial state is given by 
eq.(^. The results are shown at time Ih for the case of 
bo = O.Olao = O.OIM with the diffusion coefficients given 
by D = 10~^m^ /s and D — O.ID. The concentration is mea- 
sured in units of O.OIM. 




D [10 ' m^/s] 




FIG. 4. The diffusion coefficient of the reaction front (Df) 
and the concentration of the reaction product (co). D is the 
diffusion coefficient of the A'^ ions while D for the other ions 
is 10^'m'^/s. The units of Df and Co are 10~^m^/s and 
O.OIM, respectively. The initial concentrations are given by 
bo = O.Olao = O.OIM. The points indicated by the circles cor- 
respond to the case when the background ions have no effect 
on the dynamics. 



Having established the same scaling properties of the 
front ( pl]j2^ ) as in the case of neutral reagents, we turn 
now to the actual values of the parameters Df and cq. 
Since the motion of the reagents is modified by the elec- 
tric field (^), one expects that Df and co will depend 
not only on the properties of the reagents but on the 
properties of the background ions, as well. 

We studied the effect of the background ions by chang- 
ing the diffusion coefficient D (||) and keeping all the 
other parameters (ap, bg, D) fixed. The numerical results 
for Df and co as functions of D, are shown in Fig.^. As 
one can see, co is does not change significantly in the 
physically relevant range of 0.1 < i) < 10 (Fig.|b). The 
reason for this insensitivity of cq is that the density of 
the reaction product for ag 3> 6o and Da ~ Di, is mainly 
determined by the concentration PQl ■ 



The parameter Df is much more sensitive to the mobil- 
ity of the counter ions as shown in Fig.^a. Although the 
motion of the front is determined by the interplay of all 
four types of ions and the process is rather complex, the 
result in Figj^a can be easily understood. For ao '3> bo, 
the main effect comes from the counter ions A'^ slowing 
down or speeding up of the motion of the A^-s. If the 
diffusion coefficient D is smaller than D, the A~ ions 
are pulled back by the A~^-s (otherwise the slow A~^ ions 
would form positive charge density in the left region) thus 
fewer A~ particles enter the front which yields a smaller 
value of Df. Similar argument leads to the opposite ef- 
fect for the case oi D > D. The case {D ~ D) is special 
in the sense that the electric field (^^ vanishes and the 
result corresponds to the case of neutral reagents. 

In the next section we turn to the theory of Liesegang 
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phenomena in order to demostrate the relevance of the 
above results in the description of a relatively simple 
pattern-forming process. 



V. IMPLICATIONS FOR LIESEGANG THEORIES 

Liesegang patterns described in Sec.|l| have been much 
investigated for about a century The gross fea- 

tures of normal patterns in reproducible experiments are 
rather simple, namely the distance between consecutive 
bands Xn+i — x„ increases with band order n and the 
positions of the bands obey a spacing law p8|,pjj[: 



= 1 



Pn > 1 



P 



(24) 



where 1 -I- p is called the spacing coefScient and p > 0. 

Currently, the Liesegang phenomenon is mainly stud- 
ied as a nontrivial example of pattern formation in the 
wake of a moving front 2^ and the theories of normal 
patterns revolve around the calculation of p. The main 
feature of these theories is that the precipitate appears 
as the system goes through some nucleation- P,pO|-p4[, 
spinodal- [^, or coagulation thresholds. Most of 
these theories are rather complicated, however, and have 
been developed only recently to the level jl^jl^ that p 
can be investigated in detail and, in particular, its de- 
pendence on the initial concentrations ao and bo can be 
determined, and connection can be made to the experi- 
mentally established Matalon-Packter law . 

None of the above theories address the question of how 
the Liesegang patterns are affected by the presence of 
background ions although the existence of such an effect 
is expected. Indeed, let us take, for example, the expres- 
sion for p obtained in a simple version of the nucleation 
and growth theory [see eq.(25) in |p^] 



DrC* 



Dfico-c*) 



(25) 



where Dc is the diffusion coefficient of the C particles 
while c* is the threshold concentration of C-s. The mean- 
ing of D/ (diffusion coefficient of the C-s) and Cq (the 
concentration of C-s left behind the front) is the same 
as defined in this paper. As one can see from (p5|), the 
spacing coefficient depends both on Df and cq. Thus, 
on the basis of our results (see Fig^), we expect p to be 
affected by the background ions. 

In order to put our expectation on a firmer basis, we 
calculated p numerically employing a recent theory | [Tl| ] 
where the addition of the background ions is straight- 
forward. The main ingredients in this theory are i) a 
moving reaction front that leaves behind the particles C, 
and ii) a Cahn-Hilliard type phase-separation dynamics 
for the C particles. This theory yields the spacing law, 
and the results for p are in agreement with the Matalon- 
Packter law. Thus it appears to be a good candidate for 



the description of the Liesegang process. Since the reac- 
tion front enters the description only as a source in the 
Cahn-Hilliard equation [see eq.(3) in [|llj], one can study 
the effect of background ions by modifying the source 
according to what has been described in Sec.lV. The 



results of our numerical work for a particular case with 
bo/o-o = 0-01 (the parameters in the Cahn-Hilliard equa- 
tion were set to unity) is displayed in FigJ| 




D [lo'm^/s] 



FIG. 5. Spacing coefficient as a function of the diffusion 
coefRcient (D) of the background ions {A'^) for a case with 
bo/ao = 0.01. The circle corresponds to equal diffusion coef- 
ficients where the description in terms of neutral reagents is 
valid. 



As can be seen from Fig.^, p does depend on D and, 
actually, p can change by a factor five compared to the 
neutral case {D D) provided D decreases by a factor 
ten. One can also observe that the ionic effect is larger 
when the counterion is slower than A~ . These ob- 
sevations and the overall picture is in agreement with 
the result ( |25| ) obtained in the nucleation and growth 
theory. Indeed, cq is weakly dependent on D thus the 
main effect comes from Df. As Fig.^ shows, £)/ is a 
smooth, monotonically increasing function of D and this 
translates through eq.(p5[) into a monotonically decreas- 
ing p(i)). 

We have thus shown that the backgroung ions cannot 
be neglected in the description of the Liesegang phenom- 
ena unless the diffusivities of the ions are roughly equal. 
Although this conclusion appears to complicate the de- 
scription significantly, the reassuring aspect of the result 
is that all the complications can be absorbed into the 
parameters (D/ and cq) of the front. As a consequence, 
previous ideas about the pattern formation remain intact 
apart from the need of taking account of the renormal- 
ization of the parameters D j and cq . 
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VI. FINAL REMARKS 

A general conclusion we can draw from the present 
work is that the dynamics of reaction fronts is strongly 
altered if the diffusivities of the reacting ions differ signif- 
icantly from those of the background ions. This conclu- 
sion is based on the nontrivial density profiles found in 
a study of the the simplest reaction scheme A + B ^ C 
and assuming negligible screening length (electroneutral- 
ity approximation). We believe, however, that some as- 
pects of our results (the reaction front can still be char- 
acterized by effective diffiision constant and it still leaves 
behind a constant density of reaction product) are robust 
since they appear to follow from more general consider- 
ations and thus they should be applicable to the more 
complicated cases. 
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